Control of cardiac alternans in a mapping model with memory 
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Abstract 

We investigate a map-based model of paced cardiac muscle in the presence of closed-loop feedback 
control. The model relates the duration of an action potential to the preceding diastolic interval 
as well as the preceding action potential duration and thus has some amount of 'memory.' We find 
that the domain of control depends on this memory, independently of the specific functional form 
of the map. The memory-dependent domain of control can encompass large feedback gains, thus 
providing the first possible explanation of the recent experimental results of Hall et al. [Phys. Rev. 
Lett, 88, 198102 (2002)] on controlling alternans in small pieces of rapidly-paced cardiac muscle. 
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Sudden cardiac death, primarily caused by ventricular arrhythmias, is a major public 
health problem: it is one of the leading causes of mortality in the western world. A pos- 
sible precursor of some arrhythmias is the beat-to-beat variation of the cardiac electrical 
excitation occurring at fast heart rates . This variation appears as a sequence of long- 
short-long-short cycles (termed alternans) in important physiological characteristics such as 
action potential duration (APD) and conduction time. A focus of recent theoretical and ex- 
perimental studies is to investigate the mechanisms causing alternans and to terminate this 
response pattern using closed-loop feedback methods developed by the nonlinear dynamics 
community. Suppressing such alternans may then prevent the onset of fibrillation 2j. 

Over the last few years, several studies have demonstrated that alternans can be sup- 
p r es 8 ed with dynam ,c feedback control of the pacing Interva! QflflflQ. Control of a, t e, 
nans in the conduction time across the atrioventricular (AV) node has been demonstrated 
in both in vitro rabbit hearts \^ and in vivo human hearts The observed AV- nodal 
alternans are known to be well described by a one-dimensional map-based mathematical 
model 0. The model can be used to predict the range of control parameters that stabilize 
the desired response patterns. 

Recently, Hall et al. |4| demonstrated successful control of alternans in small pieces of 
in vitro paced bullfrog ventricles. Understanding how to control ventricular myocardium is 
important because it is the primary substrate for fibrillation. Controlling alternans in ven- 
tricular myocardium was expected to be more difficult than controlling AV-nodal alternans 
because past research suggests that its dynamics is more complicated. On the contrary, the 
experiments described in Ref. Q] demonstrated that alternans could be suppressed over a 
wide range of control parameters and over the entire range of pacing rates for which alter- 
nans was observed. The observations were compared to the predictions of two map-based 
mathematical models. Fitting the bifurcation diagrams of the mathematical models to the 
experimental data did not produce good fits for their domains of control. Specifically, con- 
trol of alternans was observed for feedback gains as large as four in the experiments, whereas 
the models predicted that the gain must be less than ~ 0.4 and limited to a small region of 
pacing rates near the bifurcation to alternans. 

The primary purpose of this Letter is to analyze a map-based model of paced cardiac 
muscle in the presence of closed-loop feedback control. The model contains some degree of 
'memory' so its dynamics displays higher-dimensional behavior. We find that the domain of 
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control can encompass large feedback gains, qualitatively consistent with the experiments 
of Hall et al. Our results demonstrate that higher-dimensional effects can enhance the 
effectiveness of control under the appropriate conditions. This may lead eventually to the 
development of methods for in vivo control of whole-heart functions. Future experiments are 
needed to determine whether the specific form of memory considered here is displayed by 
cardiac tissue. Nevertheless, our results indicate the importance of memory in the control 
of alternans. 

To set the stage for understanding the behavior of the higher-dimensional model, we first 
consider a simpler one- dimensional map that contains no memory. This mapping model is 



given by 



a simpl 



A n+1 = /(£>„), (1) 

where A n+ i is the APD generated by the (n + l) st stimulus, and D n is the n th diastolic 
interval (DI), i.e. the interval during which the tissue recovers to its resting state after 
the end of the previous (n ) action potential. We illustrate a typical bifurcation diagram 
in Fig. la. Note that the tissue has a stable 1:1 response pattern (every stimulus elicits 
an action potential of equal duration) for long pacing intervals (slow pacing rate). As the 
pacing interval decreases (faster pacing), the 1:1 response becomes unstable and a transition 
to alternans (2:2 response) occurs. At even faster pacing rates, the 2:2 response pattern 
becomes unstable and is replaced by a stable 2:1 response where an action potential is 
produced by every second stimulus. For the 1:1 and 2:2 responses considered herein, the 
pacing relation between APD and DI is D n = B — A n , where B is the basic pacing interval. 

For this simple model, the transition from 1:1 to 2:2 behavior shown in the Fig. la can be 
understood by investigating the restitution properties of the cardiac membrane. Specifically, 
to predict the pacing rates at which the 1:1 response is stable, one constructs the restitution 
curve (RC) by plotting APD as a function / of the preceding DI, as in (JTJ). Nolasco and 
Dahlen [^] proposed that the transition to alternans occurs when the slope of the RC is equal 
to unity. However, recent studies have shown that the slope of the RC at 

the onset of alternans can be significantly larger than unity and thus this criterion fails to 
predict the existence of alternans in some cases. Later in the paper, we analyze the map (JJ), 
independently of the specific form of the function /, in the presence of closed-loop feedback 
control and discuss why it fails to explain the experimental observations. In contrast to the 
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FIG. 1: Bifurcation diagrams showing existence of alternans for a) the memoryless mapping 
model (1) and b) the mapping model (2) with memory. Arrows indicate the points where slopes 
are determined. Dashed lines show unstable 1:1 response. For a), f(D n ) = A\ — A2 exp(— D n /r), 
where A\ = 392 ms, A2 = 525.3 ms, r = 40 ms [9]. For the mapping model (2), a specific form of 
the function F is taken from Refs. [14, 18], where all parameters have their typical values except 
50 ms. 
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memoryless model (JIJ, we find that a cardiac mapping model of the form 

A n+1 = F(A n , D n ) 



(2) 



can possess a domain of control that encompasses large control gains, independently of 
the specific form of F. This model has some amount of memory because it relates the 
duration of the next action potential both to the previous DI and to the previous APD, 
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tne previoi 



indicate 



i.e., the history of the dynamical system. Several previous studies 
that higher- dimensional behavior (so-called memory effects) is present in real cardiac tissue 
and thus has to be taken into account in order to correctly predict the onset of alternans. 

n 

The general form of this model was first introduced by Otani and Gilmour |15( to explain 
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empirical observations from paced dog cardiac Purkinje fibers. More recently, a specific form 
of F was derived analytically from a three-current ionic model 

Our analysis of the model (J2J) shows that it displays rate-dependent restitution so that 
there exist two primary types of RCs (the dynamic and S1-S2 RCs), which can be measured 
independently using different pacing protocols Q|. In the dynamic pacing protocol, the 
pacing interval B\ is held fixed until the tissue reaches equilibrium, and then progressively 
shortened. This yields pairs of steady-state values (A*,D*) for each B\. In the S1-S2 
pacing protocol, a premature stimulus ( "S2" ) is delivered at an interval B 2 after pacing the 
tissue with a sufficiently large number of "SI" stimuli at a pacing interval B\ so that the 
tissue reaches equilibrium. The S1-S2 RC is determined by measuring the resulting APD for 
various coupling intervals B 2 . Experimental studies have shown that the S1-S2 and dynamic 
RCs differ significantly, and have different slopes (denoted as Su and Sd yn , respectively). 
This is consistent with the predictions of the mapping model (J2J). The transition to alternans 
is governed by the combination of the slopes Sd yn and S12, so that alternans can exist when 



\F'\ 



1- (l + ^-)^2 



> 1, (3) 



where F' is the full derivative of F, with respect to A n , evaluated at a fixed point Q]. Note 
that Sdyn — 5*12 when there is no memory in the model, and that they differ substantially 
when there is large memory. A typical bifurcation diagram showing alternans in the mapping 
model (0) is presented in Fig. lb, where the slopes of the RCs are different and do not equal 
unity at the onset of alternans. 

We now consider control of alternans in both models (HI) and ©• To suppress alternans 

n 

and stabilize the 1:1 pattern, Hall et al. |4J adjusted the pacing period by an amount given 
by 

e n = -7(A n _i - A n - 2 ) (4) 

where 7 is the feedback gain. Control is initiated by adjusting the basic pacing interval B 
by e n . Applying the control technique to the mapping model (JTJ) yields 

A n+1 = f(D n ) = f(B + e n - A n ). (5) 

The linearization of (jHJ) in a neighborhood of the fixed point (when A n = A* and e n = 0) is 

e*> (6) 



A n +\ — A* + 



dAr. 



{A n ~ A*) • <)J 



de 
f-p- 



f-p- 



where f.p. denotes evaluation at the fixed point. Since the controlled pacing relation is 
D n = B + e n — A n , the derivatives in © are 



df 



OA, 



We rewrite expressions 



f-P- 
and 



df 
de 



df 



f-P- 



dD r 



(7) 



f-P- 



e n+l 

\ a n+ i J 



e-, 

V «n J 



(8) 



in matrix form as 

f * \ ( >i \( 5 n \ 

1 1 

V i o o) 

£ n /j, and a n = <5 n -i- Control is successful when the eigenvalues 
of (JHJ) fall within the unit circle in the complex plane. This occurs when \i and 7 lie within 
the region defined by the curves 
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where 5„ = A r 



A*, e n 



7 



1 + fi 



7 



(l-/i) 2 + 4 



(9) 



2// ' 2/i 
as shown in Fig. 2a. Condition © defines the domain of control for the mapping model ((TJ). 
Note that the domain does not depend on the specific functional form of /, but only on the 
value of its derivative \x evaluated at the fixed point. Alternans may exist in the uncontrolled 
mapping model (0) when jj, < — 1. The domain of control indicates the values of the gain 
7 that should be used to establish control of alternans for different values of the Floquet 
multiplier fi, which describes the stability of the map in the absence of control. Note that 
the gain necessary to establish control in the region where alternans exist in the absence of 
control is small (0 < 7 < \/2 — 1), in contradiction with the experimental observations ^| 
where 7 can be as large as 4. Very different behavior is found for the mapping model (J2J 
with memory. In the presence of control, Eq. (J2J) becomes 



A n+1 = F(A n , D n ) = F(A n , B + e n -A n ). 

Linearizing in a neighborhood of the fixed point, we find that 

dF 



A n+1 = A* + F' (A n - A*) + 



de r 



(10) 



(11) 



f-P- 



where 



F' 



dF 

0A~n 



f-P- 



dF 



f-p- 



OF 

de r . 



f-P- 



dF 



S 



12- 



(12) 
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FIG. 2: Domains of control of alternans (shaded region) for a) the memoryless model (1) and b) 
the model with memory (2). Alternans may exist in the absence of control when /i < — 1. 

As in Eq. (JZJ), we denote fi as a full derivative of F, evaluated at the fixed point. As seen 
from Eq. (J3J), the condition for the existence of alternans in the absence of control is given 
by \i < — 1. Note that the value of S\2 does not predict the existence of alternans; instead, 
it characterizes the immediate response of the tissue to small perturbations from the fixed 
point, as was shown in Ref. Hence, we expect that Syi is more relevant for determining 
the effects of control, whereas \i dictates the stability of the steady-state response of the 
tissue. 

We can rewrite expressions (@J) and (jllj) in matrix form, using the same local variables 
as in (IHI), 



/ 8 n+l ^ 

e n+l 
\ «n+l / 



fi 76*12 

-1 1 
\ 1 0; 



(13) 



The only difference between the matrices describing the controlled dynamics for the maps 
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with and without memory is the term containing the control gain 7. In the memoryless 
mapping model ([!}, the magnitude of perturbation sensitivity is proportional to the full 
derivative /1, whereas it is proportional to the slope S12 of the S1-S2 RC in the mapping 
model with memory (0). When there is no memory, S12 = —fi, whereas they can differ 
substantially when memory is present. 

By analyzing ()13|). we find that control is successful when /z and 7S12 lie within the region 
defined by the curves 



Similar to the previous model, the domain of control (jl4|) does not depend on the specific 
form of the function F. We only need to know the values of the full derivative /i and the 
slope of the S1-S2 RC Si 2 (evaluated at the fixed point) to determine whether it is possible to 
establish control. Figure 2b depicts the domain of control for the one-dimensional mapping 
model with memory © according to conditions ()14J) . Alternans may exist in the uncontrolled 
system when /1 < —1. It can be seen from Fig. 2b that the value of 7 Si 2 necessary to establish 
control for this region is relatively small (0 < 7S12 < a/2 — 1). However, the actual control 
gain 7 can be larger or smaller than this range depending on the value of S12. Our predictions 
are consistent with our earlier expectations that S12 influences the effects of control since it 
quantifies the response of the tissue (in equilibrium) to a sudden perturbation, which the 
controller attempts to cancel. 

Comparing our predictions to experiments is not possible at this time because no exper- 
iments have measured the slope of the S1-S2 RC evaluated at the fixed point l20||. Most of 

fjfl 

the experiments (see, for instance, Refs. 17, |21|) used S1-S2 pacing protocol in which the 
"SI" interval was either fixed or set to only a few different values. Instead, to make a direct 
comparison of the domain of control obtained using mapping model (j2J) to the one obtained 
experimentally in Ref. j^J, we need to know the slope of the S1-S2 RC S12 at each point 
on the dynamic RC, i.e., for different SI. The value S12 can be greater or less than one, 
depending on the specific form of F or the specific type of tissue. For example, S12 > 1 for 
the function F used to generate the plot shown in Fig. lb. Hence, the control gain 7 must 
be less than (v2 — 1) to be in the region where control is effective. However, preliminary, 
experiments with bullfrog cardiac muscle indicate that S12 can be relatively small (less than 



1. 



7S12 



(M+1)A 




(14) 



S 



0.4 and as small as 0.05) at the onset of alternans |22j and thus the control gain could be 
large. 

Thus, the experimentally measured domain of control may be consistent with the predic- 
tions of the controlled map with memory (J2J) if S\2 is truly less than one in bullfrog. The 
map without memory does not agree with experiments as demonstrated by Fig. 2a and 
as noted previously in Ref. ^J. Our analysis is the first to suggest that memory effects 
may substantially enlarge the domain of control, regulated by Si2- Future experiments are 
needed to clarify the effects of memory on control. 

A limitation of our analysis is that the extent of cardiac memory is only one previous beat 
in the model (J2J), so that the memory is short-term. However, some studies indicate that 
long-term memory effects are present in real cardiac tissue |2l| that might be described using 
a long-term memory mapping model 12, n|. In this model, an additional memory variable 
is introduced, which is assumed to accumulate during the APD and dissipate during the 
DI. Analysis of this model [12] confirms that the slope of the dynamic RC does not predict 
the transition to alternans, but it was shown in Ref. jj] that a specific form of such a 
memory model cannot explain the experiments on control of alternans. We believe that a 
generalization of the mapping model (J2J) that includes long-term memory as in Q] may 
be a fruitful direction for future investigations. 
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